Morphology of Condon Domains Phase in Plate-Like Sample 
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Based on Shoenberg assumption of magnetic flux density dependence of diamagnetic moments 
which accounts for an instabihty of strongly correlated electron gas at the conditions of dHvA effect 
and diamagnetic phase transition (DPT) to non-uniform phase, we investigate the morphology of 
the Condon domains (CD) in plate-like sample theoretically. At one period of dHvA oscillations the 
intrinsic structure of inhomogeneous diamagnetic phase (IDP) is governed by the first order phase 
transitions between different non-uniform phases similar to the high-anisotropy magnetic systems of 
spin origin, and strongly affected by temperature, magnetic field and impurity of the sample due to 
the electron correlations. The phase diagrams of evolution of IDP with temperature and small-scale 
magnetic field in every period of dHvA oscillations are calculated. 
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I. INTRODUCTION 



Strongly correlated electron systems are subject of con- 
stant interest of physical community. The instability of 
an electron gas due to strong electron correlations at the 
conditions of dHvA effect resulting in diamagnetic phase 
transition (DPT) into inhomogeneous diamagnetic phase 
(IDP) with formation of Condon domains (CDs) [H, is 
intensively studied both theoretically and experimentally 
The realization of intrinsic structure of IDP is 
governed by the competition between long-range dipole- 
dipole interaction and short-range interaction related to 
the positive interface energy with typical magnetic length 
of Larmor radius rc [H , [lJ| . 

There are striking similarities between IDP in normal 
metals and other strongly correlated systems which un- 
dergo phase transition on temperature and magnetic field 
with formation of complex macroscopic patterns, e. g. 
the type-I superconductors and thin magnetic films with 
quality factor exceeding unity. The different technics, in- 
cluding the powder pattern and magneto-optic methods, 
successfully used for observation of intermediate state of 
type-I superconductors revealed a very rich structure [l5l | 
which in spite of all it complexity, amazingly reminds 
the variety of domain structures in thin film p^ . The 
list of the phenomena can be further continued by in- 
cluding other strongly correlated systems which driven 
by the corresponding parameters exhibit the analogous 
many-pattern behavior, e. g. Quantum Hall effect sys- 
tem, showing the charge density wave instability with 
formation of stripe and bubble phases fioj . spin-glasses 
[20| . itinerant metamagnets [2l|-j23l|. the Saffman- Taylor 
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instability of fluid-ferrofluid interface in rotating Hele- 
Shaw cells jlj and Langmuir monolayers [l^. At the 
background of these phenomena, the model of simple 
lamina periodic domain structure (PDS) of normal met- 
als d^-fl^ associated with CDs and used for explana- 
tion of electron instability seems to be rather simplified. 
Indeed, the new experiments on observation of diamag- 
netic domains by use of a set of Hall probes at the sur- 
face of the plate-like sample of silver Q reveal rather 
complicated domain structure which only at certain con- 
ditions, with tilting the sample relative to applied field, 
transforms into regular laminar structure similar to the 
observation of intermediate state of type-I superconduc- 
tor [S^l . Muon spin rotation spectroscopy (/xSR) being a 
powerful method for investigation of complex materials 
was successfully used for studying CDs The series 

of excellent experiments on investigation of diamagnetic 
instability by methods of /iSR spectroscopy jSj-iBi] prove 
the existence of CD phase in beryllium, white tin, alu- 
minum, lead and indium, but reveal also contradictions 
in attempts of quantitative analysis in the framework of 
existent theoretical considerations. 

The specific case of evolution of domain structure for 
two-dimensional electron gas in ultra-thin magnetic films 
with finite thickness L < Tc was considered in ^8|. It 
was shown that the electron system undergoes the phase 
transition into a domain state which can cross over into a 
modulated, or vortex, state with decreasing of the sample 
thickness. The consideration was restricted by the center 
of the period of dHvA oscillations. While it is known that 
IDP occupies only a part of dHvA period, it would be ex- 
tremely interesting to investigate the morphology of the 
domain patterns with the change of the magnetic field 
within one dHvA period. The influence of magnetic field 
on PDS at the conditions of diamagnetic instability was 
investigated in [2i| where the dependencies of the period 
of domain structure on magnetic field and film thick- 
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ness were calculated numerically. In particular, it was 
demonstrated by numerical calculations that the width 
of the domains with favorable direction of magnetization 
increases and the width of the domains with unfavor- 
able direction of magnetization decreases with field goes 
away from the center of dHvA period. The transforma- 
tion of PDS into modulated state near the critical point 
a =1 at the center of dHvA period was also illustrated. 
Unfortunately, the critical fields corresponding to transi- 
tions between different domain structures within dHvA 
period, including bubble lattice, separated bubbles and 
so on were not calculated and the phase diagrams were 
not constructed. 

In [T3I the equilibrium set for the system of strongly 
correlated electron gas at the conditions of dHvA effect 
was investigated in the framework of catastrophe theory. 
It was shown that in every period of dHvA oscillations the 
discontinuities of order parameter accompanied DPT is 
handled by Riemann-Hugoniot catastrophe implying the 
standard scenario for DPT, e.g. DPT is of the second or- 
der at the center of dHvA period, weakly first order in the 
nearest vicinity of this point and is of the first order at 
the rest part of the dHvA period. Thus, similar to other 
magnetic systems, e. g. the spin [isj and metamagnetic 
ones [23] , the DPT can be realized in rich background of 
nonuniform phases depending on the shape of the sam- 
ple. The studies of diamagnetic instability [13| are based 
on postulates of the phase theory when the irreversibility 
effects are ignored, and the results are applied for ellip- 
soidal uniformly magnetized phases in uniform external 
magnetic field. Thus, the shape effect accounted for by 
use of the conception of demagnetizing coefRcient n con- 
sists in expansion of the magnetic field range for IDP 
existence | x |< Xc{n) with increase of n. 

The calculations based on bifurcation theory allows to 
eliminate all possible stable equilibria of the strongly cor- 
related electron system, but the conclusion about real- 
ization of non-uniform states and evolution of IDP can 
be done only on the basis of consideration of global mini- 
mum of the free energy of a system for given configuration 
of domain patterns, including the energy of long-range 
dipole-dipole interaction, the energy of the interphase 
boundaries and the magnetic energy in applied field. The 
calculations become a rather complicated problem due 
to existence of non-local dipole-dipole interaction which 
can not be accounted for by use of demagnetizing coeffi- 
cient contrary to the studies performed in the framework 
of phase theory. Although, the relevant problems were 
solved in the physics of magnetism of spin origin [l8| . 
the general picture of formation of non-uniform phases is 
complicated by the electron correlations due to non-local 
magnetic interaction between electrons [27j . 

Using the methods developed in the physics of mag- 
netic materials, for a plate-like sample we calculate the 
phase diagrams related to complicated intrinsic structure 
of IDP which can realized through the formation of va- 
riety of different domain patterns, including modulated 
states in the nearest vicinity of critical point and regular 
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FIG. 1; Quality factor Q Eq. (jj]) defined for the plate-like 
sample of normal metal at the conditions of strong dHvA 
effect is plotted as a function of reduced amplitude of dHvA 
oscillations a (solid line). For a > 7r/2 the macrostructure 
of domain patterns are expected to be similar to the domain 
structure of thin high-anisotropy magnetic films of spin origin 
due to Q > 1. The dash line is asymptote of function Q = 
Q{a) a + 1 for a +00. 

parallel band domains which transferred into separated 
band domains, close-packed bubble lattices and isolated 
bubbles with changing the small-scale magnetic field x 
within one period of dHvA oscillations. We show that 
evolution of domain structure with temperature is differ- 
ent at the center of the period of dHvA oscillations and 
its other part within the range of existence of IDP. 

The paper is organized as follows. In Sec. |TT1 we intro- 
duce the model and relevant results of catastrophe the- 
ory which at proper choice of control variables provides 
a convenient tool for investigation of DPT. In Sec. IIIII 
we present and discuss the general conditions for realiza- 
tion of the DPT in the sample with taking into account 
the long-range dipole-dipole interaction and calculate 
the critical fields which characterized the non-uniform 
phases. We show that convectional scenario of transfor- 
mation of magnetic phases is influenced by electron cor- 
relations resulting in strong temperature and magnetic 
field dependency of critical fields. Finally, in Sec. IIVI we 
summarize our main conclusions. 



II. MODEL 

The properties of strongly correlated electron system 
at the conditions of dHvA effect in one harmonic approx- 
imation are described by the free energy functional 

G(?/; a, x) = G'"™-^) G(9™''\ (1) 

where 

G("™^) =acos(x-|-y) + iy^ (2) 
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Here, the small-scale magnetic field x = kiJ,o{H — Ha) 
is the increment of the large-scale internal magnetic 
field ^oH and external (applied) magnetic field fioHa, 
y is oscillating part of reduced magnetization, k — 
27rF/(/io-ffa)^ — 2tt/ AH, F is the fundamental frequency 
of the dHvA oscillations corresponding to the extreme 
cross-section of Fermi surface, AH is dHvA period and 
a = fj,Q Tiiax{dM/dB} is differential magnetic suscepti- 
bility [27j. In physical units x is of the order of ~1-10 
mT depending on the properties of the electron system, 
while fioH is ~1-10 T. The gradient term in Eq. ([2]) ac- 
counts for the short-range correlations on the scale of Vc 
C is coordinate. 
In the framework of bifurcation theory the diamagnetic 
instability can be handled by cusp catastrophe [13]. In 
case of thin film one can show that the bifurcation set in 
the a — X plane of control variables a and x is described 
by the following expression (see Appendix □ for details) 



a cos 



(V - 1- I Xc - y{a, Xc) |) = 1, 



(3) 



where the magnetization y = y{a,Xc) is given explic- 
itly by 2/ = a sin (xc + y) ■ Eq. Q defines the function 
a = a(xc) which is plotted in Fig. Ul^a). The bifurcation 
curve a = a{xc) divides the a—x plane onto two parts, the 
inner part is occupied by IDF, while the outer one cor- 
responds to uniform state. Offering a systematic way to 
study the diamagnetic instability, the catastrophe theory 
implies the different scenario of formation and evolution 
of domain patterns at the center of the period of dHvA 
oscillations and away of it. At the center of dHvA pe- 
riod, the stability criterion for CD formation is a >1 and 
theory predict the second order phase transition. The 
crossing of bifurcation set with lowing of temperature 
results in appearance of two-fold degenerate equilibria, 
thus, one would expect the stratification of the sample 
into the domains of equal widths and appearance of PDS 
due to the essential decrease in magnetostatic energy. 
The steady-state domain size is defined by competition 
between long-range dipolc-dipole interaction and short- 
range electron interaction on the scale of rc and calcu- 
lated by standard method based on minimization of total 
energy containing two terms, e. g. dipole-dipole energy 
energy and surface energy of separation of two domains 
[TH . Away of the center of dHvA period the transition is 
of the first order and takes place at lower temperature. 
At this condition the crossing of bifurcation set leads to 
appearance of additional local minimum in free energy. 
The energy corresponding to this new minimum is higher 
than the energy of the previous state. But, even in this 
case the formation of non-uniform phase can decrease 
the net energy due to decrease of magnetostatic energy. 
Thus, at the vicinity of phase boundary a probable new 
non-uniform phase consists of fully occupied areas with 
equal magnetization and depletion regions with different 
magnetization, e.g. separate bubbles embedded into uni- 
formly magnetized sample. With further decrease of tem- 
perature the energy difference between two local minima 
in free energy decreases, and further decrease of net en- 
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FIG. 2: (color online) Shown are (a) the free energy density 
as a function of magnetization y with two 
contributions {curves 1 and 2) according to Eq. ([2|; (b) the 
approximated free energy density Go Eq. ^ at the relevant 
range y <\yo \ where ±j/o = ±1.9 are two minima and (c) the 
free energy difference 5G — (^(""'■f' — Go. The accuracy of 
the approximation G*""'-'^' ^ Go is ~10"^. The calculations 
are made at the center of the period of dHvA oscillations with 
the value of differential magnetic susceptibility a —2. 



ergy can be achieved by organizing the separated bubbles 
into bubbled lattice and formation PDS of band domains 
with volume fraction linear dependent on magnetic field. 

In Sec. mil we show that in plate-like sample the IDF 
has complicated intrinsic structure and calculate the 
phase diagrams corresponding to evolution of IDF within 
the period of dHvA oscillations. 



III. RESULTS AND DISCUSSIONS 

Before proceeding with calculation we define the qual- 
ity factor for a plate-like sample according to 



Q 



92 G^" 
yy 



g2 Q{d^d) 

yy 



= 1 — a cos 2/0, 



(4) 



y=yo 



where G'^'^~'^^ = yo/2 is demagnetizing energy for uni- 
formly magnetized plate- like sample. Here, ?/o — yo(a) is 
the magnetization of two-fold degenerate uniform ground 
state given explicitly by the equation 



yo~a sin = 



(5) 
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at a >1. The function Q = Q{a) is plotted in Fig. [T] 
which shows that Q >1 for all values of a > 7r/2. Dif- 
ferential magnetic susceptibility a = a{^oH,T,Ti:)) is a 
function of magnetic field I-J-qH, temperature T and Din- 
gle temperature To- Therefore, it follows from Fig. [T] 
that in a wide range of temperature and magnetic field 
corresponding to the diamagnetic instability, the mag- 
netic properties of strongly correlated electron gas at the 
conditions of the diamagnetic instability are expected to 
be analogous to the conventional high-anisotropic mag- 
netic materials which are characterized by quality factor 
Q >1. 

Variation of the free energy G Eq. ([T]) with respect to 
the magnetization y at the center of dHvA period leads 
to the following differential equation 



a sm y — y ^ 
The first integral of Eq. ^ 



aridity = 



arl{dQyY = 2a cosy + + C, 



(6) 



(7) 



forms the basis for investigation of non-uniform phases in 
one-dimensional problems. By choosing a different inte- 
gration constant C one can obtain the solutions describ- 
ing separate domain walls, PDS and modulated struc- 
tures. The corresponding problems in physics of mag- 
netic materials were solved long ago (see, [l8|). We find 
it useful to present the analytical solutions of Eq. © re- 
lated to PDS for the system exhibiting the diamagnetic 
instability, the first, to see the close analogy with the 
physics of magnetic phenomena of spin origin and, the 
second, the periodic analytical solutions for IDP were 
not reported as far as we aware of. 

In a limit a — 1 — s-O"*", one can use expansion of RHS of 
Eq. dS]) in powers of h due to 6 < 6o = 1/6(0 — 1) <<1. 
Thus, Eq. ^ can be integrated resulting in periodic dis- 
tribution of magnetic induction b = bo sn((^/(5i. A:) with 
a period D—AK{k)Si, where sn is Jacobi elliptic func- 
tion, K is complete elliptic integral and k is elliptic mod- 
ulus. The parameter = rcy^ (1 -t- k'^)/{a — 1) is the 
domain wall width. At k —1, K —> 00 we arrive to 
the solution b = 60 tanh(^/i5i) describing the separate 
domain wall [1^. At k =0, we obtain the modulated 
structure b = Asin(goC) [13 1 where A = 6(a — 1) and 
go = r^^^/a^^ are the amplitude and wave vector. 

In case of a — 1 > 1, one can see that the main con- 
tribution is due to the first term in G'""*-''' . In this case 
the quality factor Q >1 (see, Fig. [1]), the sample is char- 
acterized by high anisotropy, and at the relevant range 
2/ <| yo I with a great accuracy the potential G^™*-''^ 
can be approximated by trigonometric function (Fig. [J) 
corresponding to two-fold degenerate equilibria y = ±yo 



Qiumf) _ 



-i/^osin^(f^) 
2 2 yo 



where 



Kn = 4a sin 



2 yo /-, 2 yo^ 

(1 — a cos — j > 



(8) 
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FIG. 3: (a) Domain wall width A and (b) characteristic 
length Ac are plotted as functions of reduced amplitude of 
dHvA oscillations a. The dash lines in (a) show the asymp- 
totes. In calculations of Ac (b) the value of Tc/L — was 
used. 



is a function of a and yo and analogous to the constant of 
easy-axis crystallographic anisotropy in physics of mag- 
netic materials. In Eq. ([8]) we omitted the unessential 
constant. Introducing a new variable Q — 7ry/2yo, we 
rewrite Eq. ([1]) in a form 



G 



1 



1 



i^osin^e + -A(ace)^ 



(10) 



with account of Eq. ([8|). The parameter A = a(2rcyo/7r)^ 
characterizing the the short-range interaction is analo- 
gous to the exchange constant for spin systems. This 
confirms the close analogy between high-anisotropy fer- 
romagnetic and the system exhibiting the diamagnetic 
instability. Minimization of free energy density Eq. (fTO|) 
allows to obtain the non-uniform state described by pe- 
riodic function 



y 



2yo . -ir .C 

sm sn( — 

TT 02 



(11) 



where 82 = k^J A/ Kq. If a much greater than 1, then it 
is follows from Eq. ([5|) that yo ~ tt, the parabolic term 
in Eq. ^ is negligible, and the solution has a simple 
form y = 2 sin""'^ sn(C/(52, fc) with 82 = Vck. At k =1, 
K ^ 00 we arrive to the solution y = 2 sinT^ tanh(C/(52) 
analogous to separate 180" Bloch domain wall. 

Similar to the high-anisotropy magnetic films in per- 
pendicular magnetic field a variety of different domain 
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patterns can form in a thin film of normal metal at the 
condition of diamagnetic instability a > ac where Oc de- 
pends on small-scale magnetic field x and is defined ex- 
plicitly by Eq. ©. Following [l3| and [ll| we cal- 
culated the phase diagrams in the a — x plane describing 
the intrinsic properties of IDP at one period of dHvA 
oscillations. In calculations we neglect the structure of 
the domain walls assuming that the size of the domain is 
large compared to the domain wall width, thus, charac- 
terizing the domain wall by specific wall energy 



4r, 
a = a 

TT 



and dimensionless domain wall width 

^[Q2_(a_i)2]i/2- 



(12) 



(13) 



We introduce also the commonly used in phase diagram 
calculations [l^ dimensionless characteristic length 



A. 



2Z' 



(14) 



where L is the width of the plate. The qualities a = a{a) 
Eq. ([HI), A = A(a) Eq. ^ and = Xc{a) Eq. ((Tl 
are characterized by strong dependencies on tempera- 
ture, magnetic field and purity of the sample through 
the differential magnetic susceptibility a (Fig. At 
the vicinity of point a = 1, the behavior of A and Ac 
signalizes about the presence of critical point of the sys- 
tem. Near the critical point, a 1-1-0"^, the spontaneous 
fluctuations become large, the length scale of fluctuations 
^ ~ J ~ (a— l)~i/2 has a power law of divergence, and the 
characteristic length Ac u ~ (a — 1)^/^ — s-O. There is 
no typical scale length except of the trivial lower (atomic 
distance) and upper macroscopic (the size of the system 
L) size scales. 

The standard procedure of calculation of phase dia- 
grams consists in consideration of different configurations 
of domain patterns and minimization of free energy den- 
sity containing three terms, e. g. dipole-dipole energy, 
the surface energy of the boundary states and the en- 
ergy of interaction with external magnetic field (Zeeman 
term) [l8| . In connection to this, there are two main mu- 
tually related to each other features that make the calcu- 
lations for strongly correlated electron gas at the condi- 
tions of diamagnetic instability different from analogous 
calculations in physics of magnetic materials. The first, 
the magnetic moments of oppositely magnetized uniform 
phases differ on values. And, the second, the energy of 
interaction of magnetization with magnetic field is not 
trivial, it is included explicitly in the first term of Eq. ^ 
and defines the magnetic field dependence of diamagnetic 
moments y± — y±{x) of uniformly magnetized phases 
within the period of dHvA oscillations due to 



y 



= -a^G^™-^') = asin(a;-f?;). 



(15) 



We demonstrate these differences by consideration of iso- 
lated band domains. Imbedding a uniformly magnetized 



band of width W and of magnetization y_ — y+ into a 
plate of thickness L and opposite magnetization we 
can write the total energy 

where the function 



2w;tan-i w+ -ln[w-2(l + u;-2)"' 



(17) 



is defined in [I^l in slightly different form, w ~ W/L 
is reduced band width and /i-t = ?/-|_(±a;). In deriva- 
tion of Eq. pB]) we used y-{x) = — j/+(— x) (see, also 
Fig.HJb)). The first term in Eq. ((T6)) is the self-energy of 
imbedded band, the second term describes the interac- 
tion of the band with uniformly magnetized matrix [Tsj 
and the third term corresponds to interface boundary en- 
ergy. Minimization of total energy of isolated band with 
respect to w results in 



(18) 



which defines the equilibrium band width w = w{x^ a) as 
a function of magnetic field x and differential magnetic 
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FIG. 4: (color online) (a) Phase diagrams in a — x plane for 
one period of dHvA oscillations and (b) magnetic field depen- 
dencies of the magnetization y± — y±{x) (solid lines) calcu- 
lated in the framework of phase theory for a =3. The char- 
acteristic fields xpDS,xi,s,Xso,xi,Q,Xbc and Xc are explained 
in the text. The sparse area in (b) corresponds to IDP. In 
calculations of phase diagrams the value of rc/L = 10^^ is 
used. 
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susceptibility a. Another equation corresponding to the 
saturation state is e*^*"'^ = |18,] which with the help of 
Eq. (fT8|) gives rise to 



47rAc 



= wdwg - g- 



(19) 



For a >2 the function fi± can be approximated by linear 
function fj,± = yo ± {yo - yc)x/xc, where yc = y{xc) is 
the value of magnetization at the bifurcation set (see, 
Fig. Hl^b)). In this case, instead of Eqs ([IH]) and HI]) we 
obtain 



1 — a;/xc = 7r[2 tan ^ w + 'whi{l + w ^)] 
'2.'n\c/yl = ln(l + w^^) + ln(l + w^), 



(20) 



where Xc = (1 - ydvoY 



is a saturation field 



which in a limit a >>1 coincides with the critical field Xc 
defined by bifurcation set ([3]). Eqs. ((20|) define a para- 
metric function a — a{xso) with w as a parameter. This 
function is plotted in Fig. 2] together with the family of 
phase boundary curves corresponding to other possible 
pattern configurations, e. g. periodic parallel band struc- 
ture, bubble lattice and isolated bubbles, calculated in a 
similar way. We adopted the indexes of characteristic 
fields [l^, e. g. Xbc is the bubble collapse field, xw is the 
bubble lattice saturation field, Xso is the saturation field 
for isolated non-interacting bands and Xhs is the strip- 
out characteristic field. The magnetic field xpos is the 
characteristic field below which the PDS of parallel bands 
exists. At the center of dHvA oscillations, x =0, the CDs 
of equal width exist. With increasing of the small-scale 
field, a;, the volume fraction of up-domains increases in 
the expense of down-domains. Close to the boundary 
X ^ a;p£)5 — 0+, the distortions of regular CD phase, e. 
g. bending, are expected. When x = + 0+, this 

structure transforms into the hexagonal bubble lattice 
which is more favored energetically. The further increase 
of the field x results in the decrease of the bubble density 
till the transformation of the lattice into the separated 
bubbles which exist as a minority phase in the range 
Xhs < X < Xbc- It follows from our analysis that the 
existence of different inhomogeneous phases is strongly 
affected by large-scale magnetic field fioH, temperature 
T and Dingle temperature Td due to the corresponding 
dependencies of the differential magnetic susceptibility 
a = a{iJLoH,T,T£i) which can be used for experimental 
studies of predicted effects. 



IV. CONCLUSIONS 

We investigated theoretically the phase characteristics 
of the IDF caused by the arising instability of strongly 
correlated electron gas in high magnetic field and low 
temperature. We show that the formation of IDP can re- 
alized through a variety of different diamagnetic phases, 
including the plane-parallel band structure, isolated non- 
interacting bands, hexagonal lattice of magnetic bubbles 
and separated non-interacting bubbles. 



The evolution of the domain patterns is different at 
the center of dHvA period and away of it. While at x =0 
the lowing of a temperature results in formation of mod- 
ulated structure in the nearest vicinity of critical point 
which with further decrease of temperature transforms 
into regular CD structure with the domains of equal 
width, away of the center x 7^0 the system of strongly 
correlated electron gas can underdo the series of phase 
transitions with formation of different domain patterns 
depending on differential magnetic susceptibility a. 

We think that the appearance of distinct complex 
structure of IDT has to affect the transport properties of 
electron gas under the diamagnetic instability, but this 
is beyond the scope of the article. The Condon domains 
provide an excellent system for fundamental research, 
and we hope that our studies will stimulate the further 
experiments on observation of electron instability at the 
conditions of strong dHvA effect and investigation of in- 
trinsic structure of the non-uniform diamagnetic phases. 
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APPENDIX: CUSP CATASTROPHE 



Bifurcation Set 



Let us introduce a new function 



1 + ai 



a2 



(A.l) 



which is a sufficiently smooth two-parametric function of 
(z, <y) C M"^, where a. — (ai, a2) is two-component vector 
which controls the approach to critical point. The func- 
tion Eq. (jA.ip describes the equilibrium properties of the 
system of strongly correlated electron gas with magnetic 
induction b = z BS the phase-state variable dependent 
on two adjustable, or control parameters, e. g. reduced 
amplitude of dHvA oscillations a and increment of inter- 
nal magnetic field x, connected to the components of the 
vector a — (ai, 0:2) by the relations 



(1 



a2 



(A.2) 



In mean-field theory of critical behavior of system [3ll |. 
the parameter ai controls the amount of ordering, or 
the value of order parameter (magnetic induction split- 
ting in case of DPT), while parameter a2 breaks the 
Z(2) symmetry of order parameter. We assume that 
the system has an equilibrium z = z^} with eigenvalue 
A = (j)z{zQ\OLa) = at a = ao, e. g. 



<^(z;a)-0 
d:,(j){z]a.) = 0. 



(A.3) 
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The Jacobian matrix of the system (jA.3[) 



dzao 



(A.4) 



has a maximal rank r = 2 at point (zq. ccq) due to non- 
zero value of the determinant 



det 



dz 



'^02 S 



1 



^0. 



(A.5) 



Therefore, the system of Eqs. (|A.3[) defines a smooth 
curve L in the joint space of phase-state and control 
variables {z, a). According to Eqs. (|A.3p this curve corre- 
sponds to the double degenerated points, passes through 
the point (zq, ao) belonging to the surface of equilibrium 
states (f>{z; a) = and can be parameterized by z due to 
non-zero value of the determinant Eq. (jA.5|) . Actually, 
in this case the existence of the two following smooth 
functions of z 



a = a{z), 

is provided by the implicit function theorem, 
defined by Eq. (jA.ip one can obtain 



For ( 



(A.6) 

•{z;a.) 



which is a semicubic parabola. 

From Eq. (jA.SP with use of relation (|A.2p recalling that 
for plate-like sample x = Xc — y where Xc is external 
magnetic field, we arrive to Eq. ([3])- 

In next section for completeness, we present the deriva- 
tion of the normal form for the system describe by 
Eq. (jA.ip . The real use of the normal form consists in 
scaling the original problem to the other one which is 
much simpler to solve. 

2. Normal Form Derivation 



It follows from Eqs. (|A.1|1 . (|A.3|) that there exists the 
equilibrium point of the system z = Q for which the cusp 
bifurcation conditions 



a,0(O;O) = 
5,,</)(0;0) = 



(A.IO) 



are satisfied. It is important to note that dzzz4>{^\ 0) 7^0. 
Thus, expanding the function 0(z; a) Eq. (jA.ip in a Tay- 
lor series with respect to z at z =0 and performing a 
linear scaling of the phase-state variable z 



77 — 6(1 + ai)z. 



(A.ll) 



ai = -2sin2(z/2) 
a2 — z ~ tan z. 



(A.7) 



The standard projection operator map the curve L C 
onto the space of control variables (a) giving rise to 
bifurcation set. The straightforward calculation based on 
eliminating of z in (jA.7P allows to obtain the following 
expression for the function 012 — 02 (oi): 

1 -t- ai 



The function 02 = 0^2 (ai) Eq- (|A.8|) defines the bifurca- 
tion curve on the plane of control variables (ai , q;2) where 
a fold bifurcation occurs. When ai,Q;2 <<1, Eq. (jA.8|) 
can be simplified resulting in standard cusp bifurcation 
set 



{2aif + {"ia^f = 0, 



(A.9) 



one can arrive to the following normal form for cusp bi- 
furcation 



/(r,;/3)=/3l-/32ry + ?7^ 



(A.12) 



where the components of the vector function (3 — (3{a) 
are the new control variables defined by the following 
relations 



pi=-Q^a2{l + aif, P2^-&^ai{l + aif. (A.13) 



The analysis of the normal form Eq. (|A.12[) which is topo- 
logically equivalent to the function Eq. (jA.ip is straight- 
forward (see, e. g. [Ill)- In particular, for |ai| << 1 
it gives the bifurcation set l|A.9p in the vicinity of triple 
degenerate point. 
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